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Abstract. Quasi-steady states with a disequilibrium between the gas phase kinetic temperature 
and the condensed phase internal temperature appear to be the norm of condensation environ- 
ments in space. Such states may constitute Prigogine dissipative structures. 


1. Introduction 


The processes of formation of dust particles from a gas phase in an open space 
environment and in an idealized laboratory enclosure may differ significantly in their 
physical characteristics. The reason for this is that unlike the latter environment, 
the former is usually characterized by thermal (and hence thermodynamic) dis- 
equilibrium. In addition, most condensation environments in space are characterized 
by appreciable degrees of ionization of the gas phase. These two effects can give rise 
to interesting modifications of the equilibrium thermodynamics of condensation 
applicable to the laboratory environment. Any such modifications remain largely 
unexplored to date, although the literature on the process of condensation in space 
today is fairly extensive. In this paper we present a brief discussion of the thermal 
disequilibrium aspect of the problem in an effort to draw attention to the space 
condensation environment as an interesting arena for application and extension of 
the ideas and formalisms of nonequilibrium thermodynamics. 

An extensive discussion of the physical properties of the condensation environments 
in space has been presented elsewhere (De and Arrhenius, 1979), wherein pertinent 
literature has also been cited. lonization effects (plasma effects) that influence the 
process of condensation have been discussed in this reference, and also by Mendis 
(1979). 


2. Thermal] Disequilibrium: A Norm of Space Condensation Environments 


An isolated cloud of gas and dust in the otherwise relatively empty space (the inter- 
cloud medium) represents a thermodynamic system in which an equality — or even a 
near-equality — of the gas phase kinetic temperature 7), and the condensed phase 
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(dust grains) internal temperature T, hardly ever obtains. Instead, thermal quasi- 
steady states are possible with 7, 4 7;,, where the rate of radiative heat loss of a 
grain is balanced by the rate of heat gain from the impinging gas molecules and the 
ambient radiation field. The existence of such states is easy to appreciate in the case 
of a cloud transparent in the wavelengths (predominantly tn the infrared) in which 
the grains emit and absorb most of the energy (Lehnert, 1970; Arrhenius and De, 
1973; De, 1977a). In this case the grains radiate essentially directly into the cold 
intercloud medium, and thus assume a temperature that 1s cooler than the ambient 
gas kinetic temperature if this is well above the temperature of the intercloud radia- 
tion field. It has been suggested, however, that if the cloud were opaque in the rele- 
vant wavelengths, then an equality between 7;, and 7, would be ensured — allowing one 
to apply equilibrium thermodynamics in treating the problem of condensation in 
such a cloud (Lewis, 1974). Lewis’s own argument is, however, based on a simple 
logical fallacy. He first assumes thermal equilibrium, and eventually claims to have 
proved its existence. The logic of this argument, when properly invoked, would be 
that in the opaque interior of the cloud, the radiation field is enhanced — especially 
in the infrared wavelengths — by the trapping of radiation from the grains. This 
would cause the grains in the interior of the cloud to be much hotter than those near 
the cloud boundary, and conceivably, the grain temperature in the interior could 
approximate the gas temperature. To examine this proposition, we shall study below 
the depth-dependent temperature profile of grains in an opaque cloud. 

Figure | is a sketch showing a dust grain located at a depth / below the freely 
radiating surface of a cloud which, for the purpose of illustrative calculations to be 
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Fig. 1. Sketch of a section of a disk-shaped dust cloud with a mid-plane of symmetry. 
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presented later, is assumed to be roughly disk-shaped with a midplane of symmetry. 
The definition of a cloud surface is somewhat arbitrary, and we shall later avoid this 
arbitrariness by assuming that the gas density decreases linearly from a value N = Ng 
at the midplane to NV = 0 at the surface. The gas is assumed to be entirely molecular 
hydrogen. Outside the cloud boundary, the temperature 7) characterizes the local 
radiation field — deriving from radiation from neighboring stars (but not the cloud 
itself). This temperature is defined as that which a blackbody in the region of the 
grain would attain if it did not receive any heat from the cloud. The grains, with a 
number density n(r), are all assumed to be spherical black bodies with radius a. Under 
these conditions, we find that the grain temperature at any point in the cloud Is given 
by the convergent solution of the iteration relation: 


OTb Gary = OTS + 2CNVK(T,, — Ty 4.5) + 
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where a is the Stefan—Boltzmann constant, A the Boltzmann constant, C the thermal 
accommodation coefficient, V (= (kT,./27m)''?) the average impingement velocity of 
the gas molecules (with m the molecular mass), and the subscript 7 refers to the ith 
iteration. In the above equation the left-hand side is the heat loss rate of a grain per 
unit area of grain surface. The first term on the right-hand side is the rate of heat 
gain from the local radiation field; the second term is the rate of heat gain from the 
impacting gas molecules; and the last term is the rate of heat gain by the grain in 
question from all the other grains surrounding it. The exponential term takes care 
of the fact that the radiation from the grains at location r becomes attenuated by the 
intervening grains in reaching the grain in question. 

The method of solving Equation (1) is first to neglect the last term and obtain a 
numerical solution for 7,. Since this solution refers to a case where a given grain 
does not receive any radiation from the other grains, it gives us the grain temperature 
in a region that is fully transparent for the escape of infrared radiation. This case 
has been discussed earlier (Lehnert, 1970; Arrhenius and De, 1973; De, 1977a). 
This transparent solution constitutes 7, ,, for our iterations. We next enter 7, (1) 
in the last member of Equation (1), and solve the equation in its entirety to obtain 
Ty,(2). The iterations are then continued until the successive iteration values of Ty.) 
do not differ significantly. In this way we obtain a self-consistent solution for the 
disequilibrium steady state. 

As a numerical illustration of the above discussion, we have chosen the much 
discussed primitive solar nebula where traditionally a complete thermal equilibrium 
has been assumed ~ especially in the context of the study of formation of meteorites 
(see, e.g., review by Barshay and Lewis, 1976). The nebula is visualized as a disk- 
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Fig. 2. Gas and grain temperatures as a function of depth in the cloud of Figure 1. Also shown 
is the optical depth 7 in the cloud due to dust. 


shaped cloud around the Sun (or protosun), and we consider roughly the present 
asteroidal region (which is also often suggested to be the meteorite formation region). 
The grains are assumed to be metallic iron (with a = 0.1 uw) which, with a relative 
abundance of 107° in the nebula, 1s assumed to be fully condensed wherever the 
grain temperature is below the condensation temperature of iron. We have assumed, 
somewhat arbitrarily but roughly in keeping with popular nebular models invoked 
in meteoritic studies (op. cit.), that the temperature 7), decreases from 2000 K at the 
midplane to |/e times this value at the surface, while the density N decreases linearly 
from 10'° cm~* at the midplane to O at the surface. Other parameters are indicated 
in Figure 2. Also presented in this figure are the infrared optical depths in the nebula. 
Note that the condensation temperature of tron under the given conditions is around 
1400 K (De, 1977b), so that solid grains cannot exist near the midplane. 

While the above calculations are of a preliminary nature and several next-order 
refinements are required (some of which we have discussed elsewhere; see, e.g., De, 
1977a; De and Arrhenius, 1979), the general conclusion appears to be valid that 
even at large infrared optical depths (~ 10°—10*), thermal equilibrium is not attained. 

For the sake of completeness we should also discuss the dependence of the 7;,-T, 
relationship on the gas density. This dependence is shown in Figure 3 (calculated 
from Equation (1) without the opacity term). The densities at which a near-equality 
of 7,, and 7, obtains are rather high for known or physically realistic cosmic gas- 
dust environments — certainly for the peripheral regions of the nebula. The other 
situations in which equilibration may be approached are: (i) larger cloud dimension 
H (see Figure 1); and (ii) large gaseous infrared opacity. In the former case equilibra- 
tion will be achieved near the midplane, but disequilibrium will prevail elsewhere. 
In the second case — in the present example, for instance — the gaseous opacity would 
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Fig. 3. Degree of temperature disequilibrium as a function of gas pressure (or gas density). 


have to far exceed grain opacity for equilibration to occur near the midplane. This 
appears unlikely (see, e.g., Cameron and Pine, 1973). 

Note now that if the grains have come to exist as a result of in situ condensation 
in the nebula, then at the onset of condensation the grain temperature would be 
represented more realistically by the curve marked 7, (transparent solution) — cal- 
culated for a transparent nebula — than the curve marked T, (convergent solution), 
so that disequilibrium would be even more pronounced during early condensation. 
Furthermore, if condensation occurs as a result of radiative cooling of the gas, then 
the peripheral regions of the nebula — where disequilibrium is most pronounced — 
would be a major locale of condensation since cooling would take place more eff- 
ciently here than near the midplane. Thus, in the context of the nebular type of con- 
densation scenarios invoked in meteoritics, we may make the general conclusion 
that condensation favors disequilibrium. The fact that the actual process of condensa- 
tion of meteoritic material was a thermal cisequilibrium phenomenon seems also to 
be borne out by pertinent observations (see, e.g., Arrhenius ef a/., this volume). 

The above discussion also applies substantially to the case of condensation in 
expanding shells of supernova and nova ejecta (see, e.g., Clayton, this volume, for a 
discussion of this case). The expanding gaseous envelopes here are observed to be 
transparent in the infrared wavelengths, with a resulting temperature disequilibrium 
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of an order of magnitude (see Geisel et al., 1970, and Figure | of Arrhenius ef al., 
this volume). 


3. Condensation in Disequilibrium: Prigogine Dissipative Structures 


Dissipative structures — a concept of nonequilibrium thermodynamics introduced by 
Prigogine (see, e.g., Prigogine, 1978, and earlier references therein) — are states of 
organization of matter where disequilibrium becomes a source of order. Higher states 
of organization can arise within a physical system at the expense of energy dissipated 
to the environment of the system. As one moves from thermodynamic equilibrium, 
one expects progressive chaos; however, if one moves sufficiently from equilibrium, 
one may encounter the Prigogine states where again a tendency to order begins to 
evolve (see, also, Prigogine, 1979). 

Condensation of a vapor into dust grains represents a transition towards a higher 
state of organization or order. This is especially so when the grains exhibit a crystal- 
line morphology — as is the case for many tangible samples of space condensates 
represented by meteorites (see, e.g., Wark, 1979). We suggest here the possibility that 
condensation in space may be a dissipative structure phenomenon, having little to do 
with the equilibrium state where much of the attention of space condensation studies 
has been focused to date. 

We consider here the case of condensation onto a bulk condensed phase. We thus 
bypass, for the present purpose, the important problem of nucleation (see, e.g., 
De and Arrhenius, 1979; De, 1979; and articles by Draine and Herbst, 1979). 
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Fig. 4. Rate of condensation or evaporation of molecules as a function of the grain temperature 
in the state of thermal disequilibrium, showing the temperature domains of grain growth and the 
emergence of a Prigogine dissipative structure. 
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Progressive growth of a grain in a state of temperature disequilibrium occurs when 
the rate of influx of the condensable molecules to the grain surface times their sticking 
coefficient (i.e., the rate of condensation) exceeds the rate of efflux of the evaporating 
molecules. The former rate is determined by the gas density and the gas temperature 
while the latter rate is determined by the grain temperature, where the gas and the 
grain temperatures are coupled through Equation (1). These two rates are sketched 
in Figure 4 for a typical substance such as iron (see, e.g., Lehnert, 1970; De, 1973). 
The two curves intersect at grain temperatures 7,,, T,. and 7,3 which are charac- 
teristic of the given substance. The asymptotic increase of the rate of condensation 
near T,, is associated with asymptotic increase in the gas density, and hence an 
approach towards equality of 7, and 7;, (see Figure 3). Thus the domain of condensa- 
tion for 7, > 7,3 represents the conventional equilibrium condensation. As we move 
towards the left of the diagram, equilibrium is lost, and so is the tendency of the 
system to proceed towards a higher state of order by condensation. However, if we 
move sufficiently from equilibrium, then towards the left of 7,, we encounter a new 
state with a revived tendency to order. This is a dissipative structure. Order evolves 
(i.e., condensation takes place) as the entire dust cloud continues to lose energy to 
the intercloud medium. From our argument in the previous: section, it seems likely 
that condensation in space occurs through the emergence of such dissipative structures, 
rather than in equilibrium. 


4. Remarks 


The simple discussion in this paper suggests that a proper study of the process of 
condensation in a space environment should include any effects arising from thermal 
disequilibrium. The concept of dissipative structures may find an important applica- 
tion in this regard. Conversely, by employing this concept, interesting perspectives 
may be obtained on the condensation process and the condensation environments 
in space. 
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